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Abstract 

The integrability of a family of hamiltonian systems, describing in a particular case the motion 
of N "peakons" (special solutions of the so-called Camassa-Holm equation) is established in 
the framework of the r-matrix approach, starting from its Lax representation. In the general 
case, the r-matrix is a dynamical one and has an interesting though complicated structure. 
However, for a particular choice of the relevant parameters in the hamiltonian (the one corre- 
sponding to the pure "peakons" case), the r-matrix becomes essentially constant, and reduces 
to the one pertaining to the finite (non-periodic) Toda lattice. Intriguing consequences of 
such property are discussed and an integrable time discretisation is derived. 
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1. Introduction 

Interest in integrable Nonlinear Partial Differential Equations (NLPDEs) modelling phys- 
ical phenomena and possessing "exhotic" localized solutions has been revived in recent years 
thanks to the work of Camassa and Holm [1] . 

Namely, in [1] an integrable NLPDE, included as a special case in the general scheme 
provided by Fokas and Fuchssteiner in 1981 [2], has been rederived in the context of water 
waves, and it has been shown to possess, under suitable assumptions, a new kind of localized 
solutions, denoted as "peakons". 

On the other hand, a fruitful research activity has been in the past devoted to describe 
the time behaviour of the variables characterizing classes of special solutions of integrable 
NLPDEs in terms of integrable hamiltonian systems with finitely many degrees of freedom 
[3 - 7]. Quite naturally then, a similar approach has been pursued to investigate the N- 
"peakons" dynamics, by different authors with different perspectives, and results have been 
remarkable [1, 8, 9]. In particular, in ref.[9] it has been proven the complete integrability 
of the hamiltonian system defined, in canonical coordinates, by the following Hamilton's 
function: 

^ N N 

Ti = 77 ^ ^pjPkjX + fjco8k[u(qj -qk)] + fjf wah[u\qj - q k \]) (1.1) 

3 = 1 k=l 

which reduces, in the special case \x = (fi') 2 , to the one describing the motion of N "peakons" 
of the Camassa-Holm equation. 

In this paper, the complete integrability of the hamiltonian (1.1) will be proven in the 
framework the classical r-matrix approach [10, 11]. In this way, we will not only succeed in 
casting the system under scrutiny in a general unifying scheme, but we will also be able to 
uncover the basic features that render the pure "Peakons-Lattice" case [j? = (p') 2 (hereby 
denoted as PL) quite peculiar with respect to the general case. 

In Section 2, starting form the Lax representation constructed in [8], we derive the r- 
matrix for the general case (that will be in the following denoted as GL), which turns out 
to be of dynamical type; its elegant structure allows one to give a neat characterisation of 
the Lax pair associated with any of the flows commuting with the 7i flow. In Section 3, 
we restrict our considerations to PL, and observe first of all that its r-matrix is essentially 
constant, and is identical with the one pertaining to the finite Toda lattice; the two systems 
are thus endowed with the same Lie-Poisson structure. Guided by the above finding, we 
prove that the PL dynamics is just, in different coordinates, one of the compatible dynamics 
of the Toda lattice. Moreover, by using the factorisation (QR) algorithm [12], we construct 
an integrable time discretisation (Backlund transformation) for PL. Some possible extensions 
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are outlined in Section 4. 



2. The r-matrix for the General Lattice 

First of all. we write down explicitly the Hamilton's equations associated with (1.1): 

Pj = - o— = Pj y^Pfc{-^sinh[*/(gj - q k )] + //'sign(^ - q k )cosh.[u(qj - q k )]} (2.1a) 
aq i k 

qj = — = ^p k {\ + ficosh[is(qj - q k )} + f/sinh[z/|qr,- - q k \]} (2.16) 
°Pj k 

We shall assume both \x and fi' not equal to 0. Indeed, the behaviour of the system in 
these special cases has been already investigated by F.Calogero in a previous paper [13], and 
consequently will not be considered here. With no restriction, we can then rescale to 1 one 
of the parameters, say fi, so that PL corresponds to // = ±1. 

Following basically the notations used in [8], we write the Lax representation for (2.1) 
in the form: 

L = [L,M] (2.2a) 
where L, M are N x N matrices, whose entries read: 

Ljk = y/PjPka(qj - q k ) (2.26) 

M jk = -2Ky/pjp k a'(qj - q k ) (2.2c) 

In (2.2c) a prime denoted differentiation with respect to the argument, and the function a(x) 
is given by: 

/ \ 1 VX -1 V \ X \ fn ON 

a(x) = cosh — + psmh (2.3) 

with: 



2p 



H = -— — »; k 



1 + p 1 ' 1 + p 2 

The hamiltonian (1.1) is then the following simple function of the invariants of the matrix L: 

H= -KtrL 2 + (A + K + l)(trL) 2 (2.4) 
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As the term (trL) 2 merely generates on overall translation of the system, one can assume 
A + k + 1 = with no essential lost of generality. 

If we rewrite for a moment the Hamiltonian (1.1) in the form: 

W = ^PjPkfiqj - qk), 

with 

f(x) = X+ 1 + k[1 - a 2 (x)} 

we readily discover two important features of the dynamics, which merely follow from the 
even nature of the function f(x), namely: 

i) The N coordinates qj keep their own initial ordering (see also [1]), as we have: 

— sign(g i -q k ) = 

Hence, we can assume "natural ordering": q± < qi < ■ • • . Physically, we would say that 
the particles do not cross each other. 

ii) The conjugated momenta pj keep their own initial sign, as we have: 

— signpj = 

In particular, there are no "dynamical ambiguities" in the expression y/pjpk occurring in 
formulas (2.2). 

We further observe that Newton equations can be readily derived for the flow generated 
by (1.1). The most instructive procedure is to introduce a Lagrangian through a Legendre 
transformation: 

£(q,q)= -T<(P,q) + ^Pkqk = T-C(p,q) = ^^2qjq k G jk (q) (2.5) 

k jk 

G(q) being the inverse of the matrix (F(q))j k = f(qj — qk)- From (2.5) is natural to interpret 
the dynamics under scrutiny as a geodesic motion on a manifold characterized by the metric 
tensor G(q) [1]. 

We now turn to the evaluation of the r— matrix, which is well known to be a basic tool 
in the theory of integrable systems [10, 11]. As it has been shown for the first time in [11], 
whenever a hamiltonian system is associated with a Lax matrix whose eigenvalues - assumed 
to be simple, at least locally - are in involution for a given Poisson bracket, it enjoys an 
r-matrix representation of the form: 
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{Li,L 2 } = [ri2,Li] - [r 2 i,L 2 ] 
In (2.6) we have used the standard notations: 



(2.6) 



Li = L® 1 = ^Ue 1 <g> 1 

i 

L 2 = 1 <g> L = J^-L»l ® 

i 

where {e*}^ is a basis for the matrix Lie-algebra Q which L belongs to, and, in general, the 
coefficients will be functions on the phase space {dynamical r-matrix). 

Conversely, whenever an element r e Q ®Q exists such that formula (2.6) holds, the 
eigenvalues of L and in general its invariant functions are in involution. 
We recall that the trace form on Q allows one to identify Q with its dual (the space of linear 
functions on Q), and to consider the r— matrix as an endomorphism R on Q, rather than as an 
element of Q ® Q, so that in such a "dual" picture, eq.(2.6) induces the following Lie-Poisson 
bracket between two functions on Q: 

{f,g} L = (L,[df,dg] R ) (2.7a) 

with 

[X,Y] R =[X,R(Y)] + [R(X) J Y] (X,YeQ) (2.76) 

and: 

R(X) = ^^e^eW.I) (2.7c) 

In our case, everything is very simple. The algebra Q is just gl(N, R), so that we work 
with the standard (Weyl) basis : 

{Ejkjlm = 5jl5km 

and we write: 

{Li, L 2 } = X {Ljk, L lrn }E jk <S> Ei m (2.8a) 

jk;lm 
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ri2 = ^2 r Jk;lmEjk ® Ei m (2.86) 

jk;lm 

The basic ingredient in the evaluation of the r-matrix is the fundamental Poisson bracket 
{Ljk, Li m }, whose expression in our case reads: 

{Lj k , L lm } = -[Sjiy/pkPmiatmjajk)' + 6j my /pkPi(cxijajk)' 

-^kl^PjPm(oL mk a k j)' - 5 k m^/PjPl{oil k akj)'] (2.9) 

In formula (2.9), a prime denotes again differentiation, and we have used the shorthand 
notation ctj k = «(<?j — qk)', in the derivation of (2.9) a crucial role is played by the formula: 



(-7^ + -^)a(x)a(y) = - [(1 - (fj,') )wah[i/(x + y)] - 2n'a(x + y){signx + sigm/)] (2.10) 

Looking carefully at (2.10) we readily understand that there are two essentially different 
cases, the PL case (p 2 = 1 — > (//) 2 = 1) and the GL case (p 2 ^ 1). 

A comparison of formulas (2.6) and (2.9) yields the following expression for the r— matrix: 

ri2 = -^py^signfa- - qk)Ejk ® (E jk + E kj ) 

jk 

+ 1^ " P 2 ) E - E kj ) (E fl + E tj ), (2.11) 

jkl 

where the symmetric matrix S is determined, up to functions of L, by the linear system: 

[S,L] jk = sinh[^(^- - q k )} := A jk (2.12) 

It is instructive to consider the dual picture as well, where the endomorphism R(X) (2.7c) 
reads: 

R(X) = j(p[(X^). - (*<->)+] + 2(l-p 2 )[S,X(% (2.13) 

In (2.13) by X^ we have denoted the symmetric part of the matrix X, and the suffix +(— ) 
stands for the strictly upper (strictly lower) triangular part. To write down (2.13) we took 
advantage of the natural ordering, replacing sign(gj — q k ) by sign(j — k). 

It is now a straightforward task to write down the Lax representation for the commuting 
hamiltonian flows generated by the invariant functions J 7 ^ = j^-trL k+1 . In fact we have: 

rlT 

— = {^ fc ),L} = [L,M«] (2.14a), 
ot k 
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with: 

fe-i 

M (k) = R ( L k^ = ^[ p ((L fc )_ -(L fc ) + ) + (p 2 -l)Vi 3 iL fc - H (2.146) 
8 — ' 

i=o 

Note that a remarkable simplification in (2.14a) can be achieved thanks to the identity 

fc-i 

[L^L J AL k - 1 ^] = [L k , A] . 



The PL case and its connection with the Toda lattice 

As we have already noticed, the PL case corresponds to the special choice p 2 = 1, and 
is thus endowed with the r-matrix: 

ri2 = -^py^signfa- - qk)E jk ® (E jk + E kj ) 

jk 

or, in the dual picture: 

R PL (X) = ^(p[(X^)_-(X^) + ] 

which is just one of the r-matrices that can be associated with the finite nonperiodic Toda 
lattice [14] . Such striking coincidence strongly suggests to give a closer look to the PL case, 
to uncover possible further connections with the Toda system. 

To this aim, we first notice that the inverse of the L matrix (2.2b) is in this case easily 
computable. 

Clearly, the problem reduces to the evaluation of the inverse of the matrix 

C jk = a jk = exp(±^|^- - q k \) (3.1a). 

If we introduce the quantities: 

ej = exp(±^(o i+ i - q 3 ) (3.16) 
and take into account the natural ordering, we readily see that: 

Cj,j+r = e j e j+l--- e j+r-l (3.1c) 

so that the matrix C can be explicitly inverted yielding a tridiagonal Jacobi matrix. The 
final result is the following: 
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L = 



(B 1 Ai 

Ai B 2 A 2 ... 

'•• 

V : ... A N _i B N ) 



(3.2) 



where: 



Ai = 1 ^ 



PiPi+il-eJ ' Sj Pi (l-ej_i)(l-ej) 
with eo = eAr = 0. 

The isospectral flows of L -1 are of course characterized by the dynamical equations: 



1 1 - ej_ x ej 



(3.3) 



dL- 1 
dt k 



[L-\R PL {L k )\ 



and in particular we have: 



{H, L -1 } = [L- 1 ,^^)] 



The above equation can be trivially but usefully written in terms of the new dynamical 
variable A = L _1 , taking the form: 



A = [A, ^(A" 1 )] = [A, R PL (V A lnrfetA)] (3.4). 

Thus we have the following Proposition: 

"When the variables (q,p) evolve according to the Hamiltonian (1.1), the "Flaschka 
variables" (A,B), defined through (3.3, 3.1b) evolve according to the Hamiltonian IndetA" . 
In this sense, the PL lattice is just one of the commuting flows of the (finite, nonperiodic) 
Toda hierarchy. 

Encouraged by the above property, we have further exploited the connection with the 
Toda system to derive an integrable discretisation of the PL system, i.e. a symplectic map 
preserving the conserved quantities of the continuous-time model. Like in the Toda system, 
the basic tool is the factorisation of the Lax matrix [12, 15]. Taking care of the natural 
ordering, and assuming for the sake of simplicity the + sign in (3.1a, b) it can be written as: 



L = J2 p j E n + ^\fWPk e M^{q.k-q3)}Ejk + Yl VWPk^M^ij -ik)]E jk (3.5) 

j k>j " k<j 

It is then possible to find out explicitly (and uniquely, up to signs) an upper triangular 
matrix U(L): 
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U(L) = ^2^2u jk E jk (3.6a) 

j k>j 



such that: 



L = U{L) x [U(L)Y (3.66) 
where the superscript t denotes transposition. In fact one finds: 

u )k = Pj[exp[u(q k - qj)] - ex.p[is(q k+1 - qj)] (3.6c) 

The above expression holds for 1 < k < N, provided we set qN+i = —oo. 
The new (or "Miura transformed") L, say L, is then given by: 

L = [U(L)Y x U{L) (3.7a) 
which yields for the entries uj k (3.6a) the alternative expression: 

u ]k = Pfc[exp[z/(g fc - qj)] - exp[v(q k - ^_i)] (3.76) 

Identifiying (3.6c) and (3.7b) we get, after some trivial manipulation: 



Pk _ Pj 



exp[u(q k - q k )] - exp[v(q k+1 - q k )} exp[u(qj - qj)} - exp[z/(g J - - qj-i)] 



(3.8) 



The r.h.s. and the l.h.s. of (3.8) must then be independent of the label, i.e. they have to be 
equal to the same constant, say (3. Therefore, eq.(3.8) is in fact equivalent to the canonical 
transformation: 

p k = (3(exp[u(q k - q k )} - exp[i/(g fc+1 - q k )}) 

p k = (3(exp[u(q k - q k )\ - exp[i/(g fc - q k - X )\ (3.9a) 
whose generating function is: 

S(q,q) = -^2(exp[u(q k - q k )} - exp[v(q k+1 - q k )} (3.96) 
1 k 

The canonical transformation (3.9) can be interpreted as a discrete-time flow, by setting: 



(Pj, Qj) = (pj(n),qj(n)); {pjAj) = (Pj(n + l),Qj(n + 1)) 

9 



and written in the following discrete Newton form in terms of q = q(n), q = q(n + 1), q = 
q(n-l) : 

exp[u(qj - qj)] - exp[v(q j+1 - %)] = exp[v(qj - qj)] - exp[i/(qr,- - ^_i)] (3.10) 

Moreover, as in the Toda case, such discrete flow can be explicitly integrated. Namely, if we 
denote L = L(p, q) = L n , L = L(p, q) = L n+1 , we get [15]: 

L n = W(L")L [W(L-)]- 1 (3.11a) 

From (3.11a), looking at n as at a discrete time, we can read off the interpolating hamiltonian 
T{L) defined through: 

whence it follows: 

F{L) = trLQnL-I) (3.116) 

In connection with the above result, it might be worthwhile to recall that interest in integrable 
time-discretisation has been recently revived by a number of stimulating papers [16]. 

Concluding remarks 

We would like to summarize here the main results derived in the present paper. 

First of all, we have found the proper r— matrix formulation for the family of hamiltonian 
systems recently investigated in [8], containing as a special case the "peakons" system: we 
have provided an alternative proof of complete integrability and have exhibited a Lax rep- 
resentation for the various commuting flows. Secondly, we have noticed that the PL system 
enjoys the same r— matrix as the finite nonperiodic Toda system: this property led us to 
identify the PL system with one of the commuting flows of the Toda hierarchy and paved the 
way for the construction of an integrable time-discretisation. 

Several questions are still open. The most natural one refers perhaps to the algebraic 
nature of the hamiltonian system corresponding to the special choice p = in (2.3) (yielding 
/j,' = in (1.1)), already shown to be super-integrable in [13]: namely, is the "r— matrix" 
pertaining to this special case still a "good" one? We did not check Jacobi identity in 
this special case, but anyway observe that the corresponding Lax matrix is of rank 2, and 
thus one of the conditions which guarantee complete integrability is not fulfilled (in fact, 
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complete integrability has been proven in [13] in a different way, not relying upon the Lax 
representation). A further open problem is the construction of an integrable discretisation 
for the whole family of hamiltonian systems, not only for the PL system: finding an explicit 
factorisation is however a rather hard task in the GL case, and so far we have not succeded. 
Finally, the existence of a bi-hamiltonian formulation is something to be investigated: in this 
context, the well known results available for the Toda system may be exploited in the PL 
case, while again the GL system needs a separate investigation, as the dynamical nature of 
the r— matrix makes the results derived for instance in [14] not immediately applicable. 
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